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>• ■ Abstract: We give a unified description of self-dual SU (2) gauge fields on tori of size 
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transformation, extended to the case of twisted boundary conditions. We show 
how torus calorons {It/h small) are Nahm dual to the torus instantons {It/h large). 
^ ' Holonomies are dual to the locations of constituents, this duality becoming exact in 
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1 Introduction 

The study of SU (N) gauge fields on the torus has a long history (see Ref. [||] for a recent 
account). It was initiated by G. 't Hooft [Q] who pointed out additional non-trivial topolog- 
ical features, associated to what is known as twisted boundary conditions (tbc), describing 
5^ ! the presence of center-valued conserved fluxes for non-abelian gauge fields. In particular, 
the study of self-dual gauge field configurations on the torus has become a challenging 
problem. Both within the finite volume Hamiltonian and the finite temperature formula- 
tions of non-abelian gauge theories they play an important role, and the main result of this 
paper is a rather surprising detailed dual relationship between these solutions. 

For non-zero twist, numerical methods based on the lattice formulation of gauge theories 
strongly suggested the existence of non-abelian self-dual solutions for all torus sizes |Q ||. 
Indeed, existence of smooth solutions with one unit of topological charge (Q) and any non- 
trivial twist has been established rigorously 0, even though one can prove non-existence 
for zero twist 0. (For higher topological charges existence has been proven by Taubes 
earlier). To date, however, no analytic expressions have been found for these non-abelian 
solutions on the 4-torus. This contrasts with the situation for (or S^), where one has 
the algebraic ADHM construction M. 
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An important technique for studying self-dual gauge fields is the Nahm transforma- 
tion which is particularly simple for the torus. It maps SU{N) self-dual solutions of 
charge Q on the torus to SU{Q) self-dual solutions of charge on the dual torus. The 
mapping is an involution, i.e. applying the Nahm transformation again brings one back to 
the original gauge field configuration. The only analytically known example of the Nahm 
transformation on the torus, based on a particular abelian solution |ffOl, maps it onto it- 



self [|11[]. Nevertheless, the Nahm transformation is believed to be an important ingredient 
in leading to analytical results for the self-dual solutions. 

A few strategies have been used to improve our understanding of this problem. One 
follows from the consideration of mixed compact and non-compact directions T" x 
Under the duality transformation the non-compact directions collapse and the dual space 
is an n-dimensional torus, thereby achieving a dimensional reduction. The most dramatic 
example is that of = 0, where the dual space collapses to a point, being the reason that the 
ADHM construction ^ is algebraic. The n = 1 case is relevant at finite temperature [0, |T^ 
for which the associated geometry is x S^. The case n = 3 was considered as a way to 
circumvent the no-go theorem on Q = 1, zero-twist configurations |[T^, T^|. This geometry 
is relevant for the Hamiltonian formulation of gauge fields on the spatial torus, with the 
non-compact direction identified with time. The possibility of having Q = 1 configurations 
introduces an important additional simplification, since its Nahm-dual connection is an 
abelian gauge field. 

Applying Nahm's transformation in the non-compact case demands some modifica- 
tions |TT|. In particular, there appear a finite number of singularities where the Nahm 
transformed gauge field is non-selfdual. These occur when the holonomy of the original 
self-dual gauge field (extended to U{N) by adding a curvature free abelian gauge field, 
which parametrises the coordinates of the dual torus) has a trivial eigenvalue. The holon- 
omy is defined by the Wilson loops winding non-trivially in the compact directions, when 
taken to infinity in the non-compact directions. 

An important example of the appearance of these singularities in the Nahm transforma- 
tion is the case of calorons |T2|, , instantons at finite temperature. It was already known 
from the early work of Nahm that the dual formulation involves a [/(I) gauge field on 
the circle with suitable singularities. It was only recently that explicit expressions for the 
Q = 1 self-dual configurations on x with non-trivial holonomy were obtained []T6|, IT? 



In particular it was shown how the location of the singularities in the dual field is related to 
the holonomy (the Polyakov loop at spatial infinity). These non-trivial calorons reveal the 
presence of constituent BPS monopoles, hidden at trivial holonomy due the masslessness 
of one of the monopoles. 

Another example recently analysed in detail is the Q = 1 instanton configuration for 
X R with twisted boundary conditions in time There the location of the singu- 

larities is determined by the holonomies (Polyakov loops in the three directions of T^) at 
infinite time. Away from these singularities, the abelian Nahm transformed field is self- 
dual and therefore satisfies Maxwell equations. The singularities act as sources and can be 
interpreted as point-like dyons carrying (equal) electric and magnetic charge. 

It is interesting to ask oneself how these caloron and instanton configurations are ap- 
proached from configurations on a torus. It x l^, when sending either Ig or It to infinity. 
When finite, these configurations will be referred to as torus calorons and torus instantons 
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respectively. Both types of configurations have been analysed in the past by means of 
numerical methods. In it was seen how some of the x R configurations are the 
limit of self-dual configurations on the 4-torus with twisted boundary conditions in time, 
which are known to exist |p. The Q = 1 solution tunnels between vacua where some of the 
holonomies have opposite signs. The non-existence of the Q = 1 solution without twist was 
understood as an obstruction for tunnelling between two vacua with the same holonomy. 
The situation seems quite similar to the 0(3) model on the cylinder, which turned out to 
be tractable analytically [|l^ . 

The case of non-zero spatial twist was studied numerically even earlier. The basic 
building block is a lump -the twisted instanton- carrying fractional topological charge 
{Q = ^ for SU{2)) and which is well localised in time [^, ^, |2D|. Single twisted instanton 
configurations arise for non-zero (so called non-orthogonal) twist, both in space and time, 
and they have a well-defined limit as — > cxo. Only their overall position is a free parameter, 
so that in particular their scale is fixed (to Ig). Higher topological charge configurations 



look as superpositions of these twisted instantons, distributed in the time-direction |22 



These Q = I lumps also feature prominently in numerical studies that address how 4-torus 
configurations approach the non-compact geometries of x . 

Finally, caloron configurations on a torus have been recently studied on lattices 
typically having Ig/k = 4 or higher. These torus calorons fitted extremely well to the 
infinite volume analytic results. Interestingly, twist in the time direction enforces the 
holonomy (when taking the limit Ig oo) to have zero trace, which implies the constituent 
monopoles of the caloron have equal mass. For spatial twists, however, the numerical 



analysis of Ref. |^ suggests that the constituent monopoles are forced to sit at fixed 
relative positions on the torus, but can have arbitrary masses (only constrained by the 
total topological charge). One can also study one constituent monopole in isolation by 
having suitable twists, both in space and time, such that it supports the Q = ^ twisted 
instanton mentioned before. Indeed in the limit It/lg —>■ this solution becomes time 
independent and thus has to approach a single BPS monopole. Hence also in this limit the 
twisted instanton - as the BPS monopole - is the basic building block, but now arbitrarily 
placed in space. 

Quite remarkably, as will be discussed in this paper, torus calorons are dual to torus 
instantons, with either their space (for calorons) or time (for instantons) positions deter- 
mined by the holonomies of their Nahm dual. In the limit It —>■ oo, relevant for x R, the 
Nahm dual has the limit /j — > (setting Ig = 1 = l/h for convenience). This means that 
the non-abelian cores of the torus caloron shrink to zero and one is left with the abelian 
field discussed before. Quite in the same way one can interpret the abelian dual of the 
caloron in the limit Ig oo (here setting units by It = 1, the Nahm dual has Ig —>■ 0) as 
the limit in which the instanton cores have shrunk to zero. In both cases the non-abelian 
field is self-dual at all values of k and Ig, and the violation of the abelian self-duality in 
the limiting case is thus due to ignoring the singular non-abelian cores of the constituents. 
Such understanding of the singularity may prove important for future analytic work in the 

X R case. 

In our analysis, two ingredients play an important part. One is that for lattice 
techniques have been developed p5| that allow us to study the Nahm transformation 



numerically. This was tested in a non-trivial example for SU{2) with charge 2, which maps 
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to itself. This demonstrated that the method is precise enough to become a useful tool in 
numerically analysing the properties of the Nahm dualities. This charge 2 configuration 
was actually composed from four Q = ^ twisted instantons. As we will see, it is these 
twisted instantons that are mapped on to themselves as they only have their position as a 
free parameter. 

The other ingredient is the recent formulation of the Nahm transformation in the pres- 
ence of twisted boundary conditions The Nahm transform of an SU{N), charge Q 
self-dual configuration is an SU{NoQ) configuration with topological charge Q = N/No, 
where Nq is an integer depending on the twist. If the original twist is non-trivial the Nahm 
transform lives on a torus with non trivial twist. In essence the construction is related to 
recognising twist as (for SU{2)) half-periods for a larger torus, to which the usual Nahm 
transformation can be applied. Remarkably, the dual configuration admits "half-periods" 
as well, and the dual gauge field can thus be formulated in terms of gauge fields with 
twisted boundary conditions. It is this that allows us to demonstrate that the Q = | 
SU{2) twisted instanton is mapped to itself under the Nahm transformation. 

Alternative to the formulation in terms of "doubling", one can introduce a suitable 
flavour group to compensate for the twist [^, with which the Nahm transformation can 



then be performed. Again one might study half-periods in the resulting dual torus and 
reach a self-dual configuration in a smaller torus (the Nahm-dual torus) and with non-zero 
twist (the Nahm-dual twist) Furthermore, it was shown that both formulations agree. 
From the analytical point of view this method is more general and allows a simple and 
compact characterisation of the Nahm-dual torus and twist. This will be presented in an 
Appendix. 

In this paper we will make a synthesis of all the results mentioned above, leading to 
additional insight in self-dual gauge fields on the torus and their moduli spaces. In Section 
2 we will set up the formalism, whereas section 3 shows how the torus calorons and torus 
instantons are related under the Nahm transformation, illustrated explicitly by numerical 
examples. We show that taking the appropriate limits, as discussed qualitatively in this 
introduction, reproduces the known analytic results. 



2 The Nahm transformation with twist 

Let us consider a 4 dimensional torus. This can be defined as -R^/A where A is a rank-four 
lattice of points. Let e'-"^ stand for four vectors {a G {0, 1, 2, 3}) which can be chosen as 
generators of A. We will take in general a hypercubic torus It x li x I2 x l^, for which 
e(o) = (/i, 0,0,0), eW = (0,/i,0,0), e^^) = (0,0,^2,0) and e^^) = (0,0,0,/3). Often we will 
take h = h = h = h- For obvious physical reasons we will refer to the direction with 
length It as time and to the other three as space. 

Now consider SU{N) self-dual gauge fields A^{x) defined on this torus. They can be 
looked at as fields defined on R^, periodic modulo gauge transformations under translations 
by vectors a in the lattice A. Given any element of A, a = e(s) = SaC^"^ (sa E Z), we have 
(our convention throughout this paper is that A^{x) is hermitian) 

A^{x + a) = [na]A^{x) = na{x)A^{x)nl{x) + ina{x)d^nl{x) , (1) 
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where ^la{x) form a family of SU{N) matrices, known as twist matrices. Since the lattice 
A is abelian, the matrices ^la{x) must satisfy the following consistency conditions 

ilaix + b) ntix) = e^-^-^^'*-'') + a) ilaix) . (2) 

Because the matrices belong to SU{N), the factor exp{2TriJ\f{a, b)) must be an element of 
the center Zj^i. Hence, A/'(a, b) is an antisymmetric bilinear form defined by its action on 
the generators 

Ar(e("),e(/^)) = ^ , (3) 

where Uajj is an antisymmetric matrix of integers defined modulo A^, known as the twist 
tensor. We will frequently refer to its elements in the form of two 3-vectors of integers: 
ki = rioi and rrii = ^eijkUjk- 

An important role in what follows will be played by the lattice Aq, a sublattice of A. 
This is given by the elements a G A, such that the associated gauge transformation ^^(x) 
commutes (in the sense of the composition in Eq. (^) with all other ^^(x). Equivalently 

Ao = {aeA\ Af{a, b) e Z, V6 G A} . (4) 

The quotient group A/Aq is a finite abelian group. Using the freedom to select a basis of A 
such that has the canonical Frobenius form, one can show that the order of this group 
is the square of some integer, N^. The integer A^^o plays an important role in what follows 
and depends only on the form M . It is clear that for trivial twist vectors k = rh = mod 
{M{a, b) e Z) one has Aq = A and A^o = 1- 

Now let us proceed to introduce Nahm's transformation |^. For that purpose add 
a 4 parameter family of curvature free abelian gauge fields to the original self-dual gauge 
field, A^{x) + 27iz^I. This family of U{N) gauge fields is still self-dual, as the curvature 
is not affected by the addition. The zero-modes of the Weyl equation in this background 
satisfy 

(D-2TTiz)'^'^{x) =0 . (5) 

In this notation we associate to each four-vector a 2 x 2 matrix C = 6*0/2x2 + iC ■ a 
(with (Tfc the Pauli matrices). The covariant derivative is given by 9^ — iA^{x). In the 
absence of twist (/c = m = 0), the Weyl spinors satisfy the boundary conditions: 

^:{x + a) = Qa{xm{x) , (6) 



and the index theorem [^] implies that the number of solutions satisfying Eqs. 
(labelled by a) is given by the topological charge Q. Then the Nahm transformed gauge 
field is given by 

if(^)=./rf^xvl/:(x)tAvl;^(^) . (7) 



Quite miraculously this is an SU (Q) self-dual gauge field with topological charge A^. Fur- 
thermore, the gauge field is defined on a dual torus R^/A, where A is the dual lattice of 
A. 

With non-trivial twist {k, m 7^ 0) the boundary conditions shown in Eq. (|]) are no 
longer valid, since '^'^{x) takes values in the fundamental representation of SU{N), which 
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transforms non-trivially under the center of the gauge group. Two strategies were developed 
in Ref. |^ to deal with this. One involves adding flavour and the other replicating the 



original torus. Both methods were shown to lead to the same Nahm transformed gauge 
field. 

For the flavour construction one introduces the U{NfN) self-dual gauge field A^{x) ® 
iNfxNf The enlargement of the colour space allows us to modify the twist matrices as 
follows: ^la{x) — > ^a{x) ® T*{a). The Nf x Nj constant matrices r(a) satisfy 

r(a) r(aO = e'"'^("'"') T{a') T{a) . (8) 

In this way the new U{NfN) matrices ^la{x) have no twist, and the standard formalism 
can be applied to them. The existence of solutions to Eq. (P), known as twist-eaters, was 
studied in the past P3|. The constant matrices r(a) are known to exist provided Nj is 



a multiple of Nq, the twist-dependent integer defined below Eq. (|^). The solutions form 
an irreducible set provided the matrices are precisely Nq x Nq. Furthermore, the solution 
is unique modulo similarity transformations and multiplication by constants |T|. If 
we restrict to matrices belonging to U{Nq) these constants are necessarily phases (a U{1) 
representation of the lattice A). Hence, we will consider Nf = Nq and make a particular 
choice of solutions r(a). In the appendix we show how different choices of the twist eaters 
correspond to translations in z of the Nahm transformed field. 

We have reduced the problem to the situation without twist, for a U{NqN) self-dual 
gauge field defined on the torus R^/A with topological charge = NqQ. One can then 
construct the Nahm transformed gauge field in the standard way. It will be a U{N) 
self-dual gauge field with topological charge NqN, living on the dual torus R^/A. The 
construction is in terms of the solutions of the Weyl equation in the background of the 
U{NqN) gauge field (adding 2t[z^I). However, given the tensor product form of this gauge 
field, these zero-modes are seen to be solutions of the original U{N) Weyl equation (|). 
But the boundary condition, Eq. (|), is replaced by 

vl/r(a: + a) = fi,(x)r^(a)vl/f(x) , (9) 

where i,j are flavour indices (taking Nq values). The index theorem tells us that the label 
a, running over the linearly independent solutions, takes N = NqQ values. Choosing these 
solutions to form an orthonormal set, we can express the dual gauge field as follows: 

Af{z)=tY.ld'x%-{x)^—%^{x) . (10) 



i=i " ^^/i 

?4 



Although this field is in principle living on the torus -R /A, in Ref. ||2^ it is shown that it 
possesses additional periodicities. Hence, one can define the Nahm transform of the original 
gauge field with non-trivial twist as the restriction of the gauge field given in Eq. ([To|) to 
the minimal torus, in which case it has non-trivial twist as well. It can be shown that 
this minimal torus is precisely R^/A, where A = Aq, the dual lattice of Aq (see Eq. (H)). 
This definition preserves the property that the Nahm transformation is an involution, and 



reduces to the standard construction for trivial twist. Proofs are given in Ref. [26]. A more 
elegant, basis independent derivation is given in the appendix. 
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Now we will explain briefly what is the essence of the alternative construction, which 
is directly related to our numerical implementation of the Nahm transform. We can view 
Afj{x) on R^/A as being defined on the 4 dimensional hypercube formed by the unit cell 
(spanned by e*-^^) of A, with appropriated twisted boundary conditions. Consider now a 
sublattice A' whose unit cell is obtained by duplicating the original hypercube in various 
directions, in such a way that the net twist vanishes. So one defines Ai^{x) on R^/A', such 
that its boundary conditions do no longer carry twist. A trivial way to achieve this is by 
taking A' = Aq, whose unit cell contains Nq unit cells of A, but this is not a minimal 
choice. 

To determine A' it is simplest to make a choice of generators for A such that n^/j 
takes the Probenius standard form with non-zero entries defined by no2 = —n2o = Qi and 
nsi = — ni3 = q2. We introduce pi and p2 as the smallest positive integers such that 
QiPi/N G Z. One duplicates the hypercube {pi — 1) times in either the or 2 direction and 
{p2 — I) times in either the 1 or 3 direction such that Nq = piP2- Further freedom arises 
in case Pi has non-trivial prime factors. The new gauge field on i?^/A' has topological 
charge NqQ and no twist. Now one can apply the Nahm transformation, which maps to a 
SU{NoQ) self-dual gauge field of topological charge N, defined on R'^/A'. 
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Figure 1: We present a cross section through the 0-2 plane of the various lattices that 
appear in the Nahm transformation with twisted boundary conditions no2 = —''^20 = 1 
on j A. Units are such that l\ = l\ = 5, for which R^/Aq is of the same size as the 
original torus. In this case Nq = 2 and the two choices for the minimal lattices with 
twist are denoted by A' and A" . The dual of the lattice Aq gives the torus on which the 
Nahm transformed gauge field lives, with the same twist. Duplications on the dual side 
are identical, except for interchanging the roles of A and Aq (and less importantly A' and 
A"). When a symbol overlaps with different cells it belongs to all of those cells as a whole. 

Each choice of A' will lead to a different set of zero-modes. This can only mean that the 
resulting dual gauge fields are related by gauge transformation. Gauge invariant quantities 
have to agree when derived from the various choices of A'. But this means that the 
smallest unit cell in terms of which these gauge invariant quantities can be reproduced, is 
the intersection of the unit cells of all A', which is precisely A = Aq, the dual of Aq (in the 
same way the unit cell of A is the intersection of the unit cells of all A'). Note that Aq is 
the smallest sublattice of A, that is also a sublattice of all possible choices of A' (the unit 
cell of Aq consists of A^o unit cells of A' and therefore of Nq unit cells of A) . With the twist 
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in the standard form given above, Aq is obtained by duplicating the hypercube that defines 
the original torus R^/A, {pi — 1) times in both the and 2 directions and {p2 — 1) times in 
both the 1 and 3 directions. In Fig. 1 we elucidate this for the 0-2 plane with gi = 1 and 

Pi = 2- . 

This shows that A can be defined as an SU{NqQ) gauge field on R^/Aq and has 

topological charge N/Nq. Nevertheless, as is obvious in particular for N/Nq not an integer, 

the unit cell of Aq cannot be without twist, i.e. A is also a gauge field with twisted boundary 

conditions. For SU{2) in the case that the original twist is already in the Frobenius 

standard form (where the situation of Fig. 1 applies in the 0-2 and/or the 1-3 planes), the 

dual twist can be read off from the minimal number of duplications required to go from Aq 

to any of the choices A'. 



3 The case of SU{2) (twisted) instantons 

Having presented the general formalism, let us now concentrate on the study of SU{2) 
gauge field configurations for spatially symmetric tori. From the previous formulas one 
can deduce that the symmetric structure of the torus is not always preserved by the Nahm 
transform. Nonetheless, the dual torus is always contained within the symmetric torus 
R^/A. Furthermore, if a given torus has Is/k ^ 1, its Nahm transform corresponds to 
Is/h ^ 1- Thus, the Nahm transform of the torus instanton configurations are the torus 
caloron configurations. The space of gauge inequivalent solutions defines the moduli space. 
Its dimension is given through an index theorem to be ANQ for SU{N) gauge fields of 
charge Q on the torus. Note that QN = QN so that the dimensionality is preserved by the 
Nahm transform. Indeed, in Ref. it is proven that the Nahm transformation induces an 
isometry between the moduli spaces. 



3.1 The twisted instanton Q = 



We start with the smallest value of the topological charge, for SU{2) this is Q = The 
dimension of the moduli space is 4, so that the set of solutions is discrete up to translations. 
The value of the topological charge is related to the twist by the formula p|, ^ 



Q = --^ , (11) 

where u is an integer. Fractional values of the topological charge are attainable for non- 
orthogonal twists, i.e. k ■ rfi ^ mod A^. For SU{2) we choose k = m = {0, 1,0), with 
k-m = 1. We deduce A^o = 4 and therefore the Nahm transformation is again an SU{2) self- 
dual gauge field with topological charge |. Furthermore, it lives on a torus with generators 
^ ^here e^'^) are the dual basis vectors e(°) = (1//*, 0, 0, 0), e^^) = (0,1//,, 0,0), 

e(2) = (0,0, l/ls,0), e^^^ = (0,0,0, l/h)- The corresponding twist is unchanged. 

Now we will explain how this works for configurations with Ig/lt ^ 1. This will serve 
to illustrate the precision of our numerical techniques. One can generate self-dual configu- 
rations numerically by employing the method described in Ref. . To minimise the errors 
due to lattice artifacts we actually make use of the improved cooling technique described in 
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Ref. [|T4|, with the parameter e = 0. One can implement twist by the method of Ref. 
We have generated self-dual configurations with the twist mentioned above on lattices of 
size 4 X 16^ and 32 x 8^. Notice that in the former case Is/h = 4 and in the latter case 
|. They correspond to Nahm-dual tori for this twist. In both cases the action density of 
the configuration has a lumpy profile with a single maximum in space. The position of the 
maximum does not coincide with a lattice point, but can be estimated by interpolation. 
Then, given the interpolated position of the maximum, one can determine the coordinates 
corresponding to each lattice point with respect to it. 
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Figure 2: Comparison between the action densities of a numerical Q = \ instanton with 
twist k = m = (0,1,0), living on a 32 x 8^ lattice, and the Nahm transformation of a 
Q = \ caloron (BPS monopole) with the same twist, living on a 4 x 16^ lattice. The 
depicted values correspond to points in the 32 x 8^ lattice, along the lattice direction 
niatt = (8,4,8); the peak of the action density of the lattice instanton being interpolated, 
in lattice coordinates, to n = (15.809,3.741,3.979,4.162) and shifted to (i,^,o, j^) where 
the maximum of the Nahm transform of the caloron sits. The continuum time period has 



been set to = 1 for the caloron and hence L 



for its Nahm transformed instanton. 



What we did next was to apply the numerical Nahm transformation introduced in 
Ref. The values of z, where the Nahm transformation was calculated, were those 

corresponding to the lattice points of the dual configuration. In this way the Nahm trans- 
formation can be compared with the lattice configuration obtained on the dual torus. Fig. 2 
shows the result of such a comparison for the action density along a line on the torus. The 
agreement is remarkable. It shows the way in which our expectations work and gives us 
confidence on the precision of our numerical methods to investigate other more complicated 
situations. Note that in fact Ref. [25| dealt with the case of four Q = ^ instantons (with 
Is/h = 1) on a 12 X 6 X 12 X 6 lattice, finding that the configuration is mapped to itself 
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(up to a shift). In terms of the Q = ^ twisted instantons this can be understood from the 
fact that only its position is a free parameter. Since the work in Ref. |]2S| the numerical 
methods have been greatly improved, details of which will be reported elsewhere. 



3.2 Temporal twist Q=l 

Let us recall what is known about Q = 1 configurations with m = and k ^ 0, for the 
cases It/ Is — >■ (torus calorons) and oo (torus instantons). These will be mapped onto 
each other by the Nahm transformation and we are interested in establishing how the 8 
dimensional moduli spaces are related. 

The torus caloron configurations {It <^ h) have been studied recently by e-cooling tech- 



niques [Q. As It/ Is ^ these configurations approach some of the non-trivial holonomy 
calorons |T^, namely those associated with u = ^. The parameter uj = \uj\ is defined 
through the holonomy that in the infinite volume is defined by the value of the Polyakov 
loop at infinity 

fit ^ 

exp(27rzu; ■ f) = lim Pexp(— i / AQ(t,x)dt) . (12) 

|a;|^oo JO 

These calorons with uj = ^ correspond to equal size constituent monopoles located at two 
different points. For finite Is/k the periodic boundary conditions in space slightly modify 
the action density profiles without changing the quahtative features. Already for Ig/k = 4 



the agreement is excellent [E4]. 



Figure 3: Zero-mode density profiles for the two zero-modes of the lattice caloron (left) 
on a 4 X 16^ lattice for k = (1, 1, 1), created with improved cooling [e = 0). The profiles 
fit well to the two zero-modes for the infinite volume analytic caloron solution (shown 
on the right at y = t = 0) with u = j and constituents at yi = (2.50,0.12,0.95) and 
y2 = (1.38, —0.24, 2.67), in units where P = k = 1 {or a = I) and the left most lattice point 
corresponding to x = ^ = 0. The plots give the added densities of the two zero-modes. 

Recently also the Weyl-Dirac zero-mode for the infinite volume caloron with non-trivial 
holonomy was determined analytically |3l[] . The zero- modes are more localised than the ac- 
tion density and the agreement between the numerically determined zero-mode for the torus 
caloron is therefore expected to be even better. This is illustrated in Fig. 3, where we com- 
pare these exact infinite volume zero-modes with those that are numerically constructed. 
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Since the zero-modes are the ingredients for performing the Nahm transformation, the 
present comparison shows the numerical accuracy that was achieved. For convenience we 
have taken for this comparison the parameters that were used in Fig. 1 of Ref . ||2^ , where 
one can find the action density profiles. In this particular case /c = (1, 1, 1) and we dupli- 
cate the initial torus in the time-direction, to reach a situation without twist. This gives 
a charge 2 configuration with two zero-modes, of which one is periodic and the other is 



anti-periodic. From Ref. ||3T|] one knows that each is localised on one of the constituents. 
The numerical procedure was to combine the two zero-modes of the doubled lattice in two 
orthonormal modes, each maximising the overlap with one of the two constituents. Using 
the large degree of localisation, we have simply added the zero-mode densities together 
to show the result in one figure, but we stress each lump corresponds to the density of 
one of the zero-modes. There are no free parameters involved in this comparison and the 
agreement for the zero-modes is indeed impressive. 

In the infinite volume the parameters of the caloron are described by the spatial posi- 
tions of the two constituents (6 parameters), one overall time position and a U{1) phase 
rotation that can be undone by a gauge transformation. For large separations of the con- 
stituents the caloron becomes static with the two constituent monopoles showing up as 
separate lumps. On the other extreme when the constituent monopoles are close to each 
other they fuse into a single lump that looks like an ordinary instanton in R^. 



Torus instanton configurations {Ig <^ k) were studied a few years ago [|1^, |T5[. As 
h/k the configurations evolve into well-defined solutions on x R. One feature 
shared by all these torus instanton solutions is that they are local in time. At large times 
the configurations have to approach a pure gauge, to guarantee the action stays finite. The 
gauge fields with zero curvature on at t = ±oo are characterised by their holonomies 
(spatial Polyakov loops) 

= ^ hm P exp(-i J^' Aj{t, x)dxj) . (13) 
In the presence of a temporal twist k, the holonomies at both ends are related 

When we consider It/l^ finite, the instanton has 8 parameters, but the holonomies can 
not be rigorously defined. In the limit It/h — > oo they would, however, be related by the 
twist factors. It is likely that in this limit the three holonomies are part of the moduli 
that describe the solution (for the infinite volume caloron, gauge zero-modes associated 
with a variation of uj are not normalisable and u is not a moduli parameter). There are 



indications that the holonomies fix the scale of the torus instanton |T^, |T8[. In addition 
there are of course the four position parameters, and presumably a phase that describes 
the colour orientation relative to the twist matrices. 

The dual gauge field for instantons on x R with temporal twist was shown to be 



a solution of the abelian Bogomolny equations with point sources [0 that have (due to 
the self-duality) equal magnetic and electric charges. From this one computes the Maxwell 
field analytically by summing over the periodic copies of these charges. Essential is that 
the location of the point charges is determined by the holonomies P^. At these specific 
points on the dual torus the Nahm transformation is modified due to boundary terms in the 
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non-compact directions. This leads to four singularities, whose locations are determined 
by requiring exp{—27riljZj)P^ to have a unit eigenvalue (this leads to a zero-energy state 
in the Weyl-Dirac Hamiltonian at t = ±00 and therefore to non-exponential decay of the 



zero-mode). The charges are all equal in absolute value, due to Dirac quantisation [jTS 



and opposite in sign for those associated to the holonomies at either end. With temporal 
twist, the location of the sources with one sign of the charge are displaced by ^(fj'j ■^) 
with respect to the location of the two sources with the opposite sign. 




Figure 4: Ratio of the action densities of the Nahm transformed gauge field for a x 8^ 
lattice, dividing the It = 40 result by that of It = 20. The (approximate) holonomies are 
given by ^TrP^" = -iTrP+ = (0.86,0.76,0.082). 

We illustrate in Fig. 4 how the singularities arise as It increases. The results are based 
on a charge 1 instanton on an 20 x 8^ and 40 x 8^ lattice with temporal twist k = (1, 1, 1) 
and holonomies defined by iTrP^^ = — ^TrP^^ = (0.86, 0.76, 0.082) (changing It and making 
sure the holonomies stay fixed can be achieved by cooling with open boundary conditions, 
introduced in Ref. [|15|). We choose units such that /<; = 1, so that the sources with positive 
charge are located at z = (0.415,0.387,0.263) and z = (0.585,0.613,0.737) and the ones 
with negative charge at z = (0.915,0.887,0.763) and z= (0.085,0.113,0.237). Fig. 4 is 
plotted in the plane defined by z^ = 0.263, which is through the position of the first source 
and still close to the position of the last source. Shown is the ratio of the action density 
for the 40 x 8^ lattice over the action density for the 20 x 8^ lattice. 

The nature of the singularities is illuminated by the observation, made before, that 
the Nahm transformation of the torus instantons are the torus calorons. Thus, the action 
density peaks correspond to the location of the BPS constituent monopoles. Their mass 
(associated to w = j) should be given by Sn'^lt (in units where Ig = 1). The result 
is compatible with the fact that the peak heights of these BPS monopoles scale as l^. 
Away from the core of these BPS monopoles the field becomes abelian and should become 
independent of It, as the charge is fixed. Fig. 4 shows that indeed the ratio becomes one 
outside of the core of the monopoles. The non-abelian core of the BPS monopole shrinks 
to a point in the limit It 00, and one is left with the resulting abelian field. This is 
illustrated in Fig. 5, where we plot the square of the electric field in the plane defined by 
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Z3 = 0.5. The numerical result is based on the 40 x 8^ lattice, for which the plane 23 = 0.5 
is far enough away from all sources not to be affected by their non-abelian dressing as 
SU{2) BPS monopoles of mass 407r^ (in units where = 1). The agreement, that involves 
no free parameters, is indeed very good. 



Figure 5: Comparison of E'^ in the plane Z3 = 0.5 for the Nahm transformation of the 
40 X 8^ instanton configuration described in the caption of Fig. 4. On the left is the 
numerical result and on the right the analytic result for x R. 

The Nahm transformation was constructed by doubling in time, such that at finite k 
this corresponds to a charge 2 torus instanton, which is mapped to a charge 2 torus caloron 
defined on the hypercube [0, x [0, 1]^. Our general formalism tells us that this is formed 
by duplicating the Nahm-dual torus, which is half the size. The latter corresponds to a 
skew lattice, since the twist is not in the Frobenius standard form. The additional lattice 
generator is precisely the vector i(0, j^, ^, 1^) giving the displacement of the 2 positive 
versus the 2 negative sources. 

In summary, we have illustrated how the Nahm transformation maps torus instantons 
into torus calorons, and how the (approximate) holonomies are mapped into the positions 
of the constituent BPS monopoles. This relation is exact in the limit that It —>■ 00 (resp. 
0). To test its validity at finite It/ls we studied the case with temporal twist k = (0, 1,0) 
(m = (0,0,0)). One cannot have both the original and Nahm-dual torus to be spatially 
symmetric, due to the asymmetric twist, but that is not important for the general con- 
clusion about the mapping we are studying. If our original torus has size k x l^, A = Aq 
has SIZG 77/~ X 7~ X 777— X Y~ . The twist vectors for the Nahm transform are the same as the 

Zti Ls ^I'S >'S 

original ones. 

By means of improved cooling {e = 0) we have generated a lattice torus caloron config- 



uration [^. The relative distance between the two constituent monopoles can be tuned 
by using different values of the parameter e On the other hand, it is possible to use 
a modified version of cooling |]15| (fixing the field to have zero colour magnetic fields at 



t = and t = It in terms of prescribed holonomies) in order to produce torus instanton 
configurations with fixed holonomies. If we choose these to be the ones determined by 
the location of the monopole constituents of the torus caloron, then the numerical Nahm 



transformation |25| can be applied to each of these configurations and has to give back 



the other configuration, up to a translation. This is illustrated in Fig. 6. The somewhat 
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larger numerical differences are in part due to the not sufficiently large ratio for It/h which, 
for the required values of the holonomies, enforces a rather narrow instanton which is not 
easily stabilised on the lattice. Results for a smaller value of It/h (not shown) confirm this 
improvement with increasing It/h- Going to even a larger ratio is however computation- 
ally too expensive. The comparison confirms that the holonomies of the time-twisted torus 
instanton are part of the moduli. 



18.75 — 



12.5 — 



6.25 — 



+ 



Latt. caloron : Ellipses 

Nahm(Latt. instanton) : 
Circles 



Figure 6: Comparison between the action densities of a numerical Q = 1 instanton with 
twist k = (0,1,0) and m = 0, living ona32x8x4x8 lattice, and the Nahm transformation 
of Q = 1 torus caloron with the same twist, living on a 4 x 32^ lattice. The constituent 
locations were used to fix the holonomies in generating the torus instanton. Depicted values 
correspond to points in the 4 x 32^ lattice, along the x direction. The continuum units are 
such that = 1 for the caloron. 



3.3 Spatial twist Q = 1 

We will conclude by discussing the case of spatial twist {k = 0), for which we take rh = 
(0, 1,0). If our original torus has size It x 1^, the Nahm dual torus for this twist has size 
X ^ (to obtain instead a spatially symmetric dual torus we should apply 



i X J- X 

It 2ls Is 2/ 

the Nahm transformation to a torus with size It x 



^Ig X Ig X ^Is)- 



As before, the dual 

gauge field has the same spatial twist as the original gauge field. The phenomenology of 
these configurations, for the large aspect ratios we are interested in, is different to the 
previous case with temporal twist. But there are also quite interesting analogies that may 
be important for the case It = Is, where from the Euclidean point of view there should be 
no essential difference between twist in space or time. 



The analysis of Ref. |^4[ shows that for Is/k large, torus calorons with spatial twist 
are given by configurations that approximate the infinite volume caloron with arbitrary 
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holonomy, parametrised by G [0, |], such that the two constituent monopoles have in 
general differing masses, IQir'^uj/lt and 87r^(l — 2uj)/lt. But now their relative positions are 
fixed and determined by the spatial twist. This contrasts with the situation for temporal 
twists for which the value of w is fixed, but the relative position is arbitrary. 

Torus instantons with twisted boundary conditions in space, and Is/h small, have been 
studied extensively in relation to the Hamiltonian formulation in a twisted box |^. The 



outcome of these studies |^ showed that the configurations can be described in terms of 
two Q = \ twisted instantons {2Q twisted instantons in the case of higher charges), that 
have both twist in space and time. The net twist in time, however, vanishes. We can view 
the Q = \ constituents of the Q = 1 instanton in very much the same way as the monopole 
constituents for the caloron. It was found that the space-time locations of these twisted 
instanton lumps can be arbitrary. Only when they are very close to each other, they merge 
into a single lump. In that region, when the size of the lump becomes small compared to 
Isi the configuration behaves like the ordinary instanton in i?^. There are two very simple 
specific examples where this constituent nature of the torus instantons with spatial twist 
is apparent. Take a. Q = \ instanton with twist m = = (0, 1, 0) and duplicate either in 
the or 2 direction to remove the time twist. The result is a Q = 1 instanton with spatial 
twist. Either the two lumps have the same spatial positions and are separated maximally 
along the time direction, or they have equal time positions and are maximally separated 
along the space direction. 

We now turn our attention to studying how the Nahm transformation relates these 
torus instantons and calorons. For the above example built from two Q = \ instantons, 
this can directly be understood in terms of the results of section 3.1, where it was shown 
that the dual of such a Q = 1 instanton is a BPS monopole, with twist k = m = (0, 1, 0). 
Duplication of the instanton in the direction is dual to duplication of the BPS monopole 
in the 2 direction, and we do get a caloron with maximal separation of the {u = \) 
constituents in the direction of the spatial twist. On the other hand, duplication of the 
instanton in the 2 direction is dual to duplication of the BPS monopole in the direction. 
This doubles the mass of the constituent monopole and corresponds to a caloron with 
trivial holonomy = or i (the other constituent is massless and therefore absent). This 
establishes that, with twist in space, u need not be \. 

To deal with a more general case, we consider the Q = 1 instanton with spatial twist 
being made out of two Q = i (m = (0, 1, 0), k arbitrary such that k-ffi = l(mod2)) with an 
arbitrary time separation, i.e. no longer obtained by duplication in the time direction. To 
study some of the properties of its Nahm transform, we consider the limit — * 0. In this 
limit all spatial dependence is removed and the two Q = \ instantons shrink to two points, 
specified by their positions in time. In one dimension curvature can be non-trivial only 
at point-like singularities, so the curvature vanishes everywhere, except at the locations 
of the Q = \ instantons, where the field strength becomes singular. The curvature free 
regions are described by a constant abelian gauge field Aj, making a jump at the positions 
of the Q = \ instantons. In this limit, for a single Q = \ instanton, the twist determines 
the jump in the constant gauge field to be AA^ = Txkj/lg (relevant for the case of space 
twist, so as to assure that Eq. (0) is satisfied. In the general case, where the values of the 
Polyakov loops are not fixed (see below), there is more freedom). The Q = 1 instanton is 



15 



therefore described by 

= 27r(e, + ifc,X[t„,i,](t)/y , (15) 

where (and Aq) are arbitrary constants (fixed for spatial twist, see later) and Xlta,h]i^) 
is the characteristic function of the interval t G [ta,tb]- This result agrees exactly with 
the Nahm transformed gauge field of the infinite volume caloron ||1^. We can identify 
\k/ls{= anp^) with the distance vector between the two constituents (in units set by 
It = 1) and u = ^{tf, — to) with its holonomy (see the second of Ref. [0, Eqs. (41) and 
(65), for the precise definitions of p and a). 




Figure 7: Comparison between the action densities of the numerical Nahm transform of a 
torus instanton and the appropriate infinite volume analytic caloron solution. The original 
field lives ona80x4x8x4 lattice with twist = (0, 0, 0), m = (0, 1, 0). Units are set 
by = 1. The temporal distance between the two Q = | constituents for the instanton 
corresponds to a caloron with u = (io/2 = 0.311. The plot is along the line connecting the 
centers of the two constituent monopoles, defined by zi = = 5 as determined by the 
spatial holonomies of the instanton. 

It may seem that for the Nahm gauge field of the caloron ta = —uj and tb = ui are 
fixed, but one should of course realise that A is still self-dual, away from the singularities, 
when t is shifted over a constant. We conclude that the holonomy of the caloron fixes the 
relative locations of the Q = \ instanton lumps under the Nahm transformation, whereas 
the holonomies of the instanton, determined in each of the two fiat regions (for /<j ^ 0), fix 
the relative locations of the constituent monopoles of the caloron. Their difference vector 
^m/ls, remarkably, agrees with maximal separation on the dual torus along the direction 
determined by the spatial twistQ. This derivation is as rigorous, in the limit — > 0, as the 

^Note that on the dual torus a displacement by ^ along directions 1 and 3 corresponds to a shift by 
a full torus period. Hence, a shift by ^k/lg is identical to ^rfi/ls, for any k with k ■ fh — l(mod2). 
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one for x [0, It], in the limit It ^ oo and has given a beautifully consistent picture of the 
dual relationship between the moduli spaces of the calorons and finite volume instantons. 

To test this picture at finite Ig we have made use of our numerical methods. We started 
with a configuration on an 80x4x8x4 lattice containing two Q = ^ torus twisted instantons 
(with m = k = (0,1,0)) separated in time by a distance do. In units where /* = 1 (so 
= 1/10), the lattice A is generated by Cq = (1, 0, 0, 0), ei = (0, IJ2, 0, 0), eg = (0, 0, h, 0) 
and 63 = (0,0,0,/,/2). The dual torus RyKo has as unit cell [0,1] x [0,l;^f. This 
configuration is constructed in several steps. First one starts with a. Q = ^ configuration 
ona40x4x8x4 lattice with twist k = rh = [0,1,0). Using a time reversal and parity 
transformation one obtains a second Q = ^ instanton with the Polyakov loops at t ~ ±oo 
interchanged. One can glue the two configurations together at arbitrary time separation 
(Iq and transform them appropriately to trivial twist matrices in time. If do ^ Is this 
procedure gives an approximate self-dual solution on the required torus and with spatial 
twist. The small deviations from self-duality induced by the procedure can be eliminated 
by cooling with e = 0. We end up with a. Q = 1 configuration with the desired properties. 
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Figure 8: The same as in Fig. 7, but for a temporal distance corresponding to do = 0.786. 



A few remarks are in order, concerning the values of the Polyakov loops for the space 
twisted instantons. Eqs. ( P^ and (|TE|) are only valid in the gauge where the appropriate 
f2a(x) are trivial. Such a gauge was assumed for the discussion in sect. 3.2. The general 
formula for a Wilson loop that closes by the shift over a period reads 

P^(x) = Pexp(-z / A^{x)dXf,)n^(t.){x) . (16) 

In the presence of a space twist the Pi are not arbitrary, but uniquely determined (up 
to some discrete transformations) by the twist. This can be seen as follows. Since the 
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configuration is exponentially localised in time, it goes to a pure gauge at ±00. Now one 
can choose a gauge in which Ai(t = 00) = 0, fixing = in Eq. (0), for which the twist 
matrices f2g(i) are constant and the holonomy is given precisely by these matrices. The fact 
that fig(i) and flg(3) anticommute (as a result of the twist), implies they can be brought to 
the form iai and ia^, respectively. On the other hand, the constant twist matrix for the 2 
direction has to commute with these Pauli matrices and therefore must be ±/. This fixes 
the Polyakov loops for the Q = 1 space twisted instanton at both t = 00 and t = — 00 to 
be jTrPi = ^TtP^ = and |TrP2 = il, depending on whether we choose the twist matrix 
in direction 2 to be ±/. The location of the caloron constituents are determined by the 
holonomies to he at z = (5, 0, 5), and the one displaced by ^m/lg from it, i.e. z = (5, 5, 5). 

We have generated two such configurations with time separations between the Q = ^ 
instanton constituents given by = 0.622 and do = 0.786. In Figs. 7 and 8 we compare 
the Nahm transformed results with the infinite volume analytic caloron solutions and the 
holonomy fixed to u = do/2. Plotted is the action density along the line joining the 
two monopoles. The comparison is very good taking into account that there are no free 
parameters. The mismatch at the peaks can be attributed to the fact that the data points 
correspond to a finite ratio It/lg = 1/10) which is zero for the analytic result. A direct 
comparison with numerical data at this fixed ratio is difficult due to its small value, and 
seems unnecessary given the agreement shown in Figs. 7-8. 

4 Conclusions 

By now our knowledge of self-dual configurations on the torus has come from various 
sources. Apart from general existence proofs and consequences of the index theorem, in- 
formation in this field has arisen from both numerical studies based on lattice gauge theories 
and analytical results obtained on manifolds having compact and non-compact directions. 
Most of the interest and information refers to self-dual configurations on spatially symmet- 
ric tori of size It x for large and small aspect ratios, known as torus instantons {Ig/k ^ 1) 
and torus calorons ^ 1). In this paper we have studied the Nahm transformation 

that relates SU{N) self-dual gauge fields of charge Q with twisted boundary conditions 
on a torus, to SU{NoQ) self-dual gauge fields of charge N/No on the dual torus. Here 
No is determined by the twist. We found that torus calorons and instantons are dual to 
each other. This provides information on their respective moduli spaces, which is fully 
consistent with the structure suggested by numerical studies. In particular, it allows us 
to understand the findings of Ref. ||2^. The most notable result is the duality observed 
between holonomies and the location of constituent structures (BPS monopoles or twisted 
instantons), as suggested by analysis of the Nahm transformation of configurations living 
in X R and R'^ x S^. We showed how the abelian Nahm connections with singularities 
obtained in this case, result from the collapse of the non-abelian core of the constituents 
into point-like singularities, which act as sources of the surviving abelian field. Thus, all 
the information fits nicely into a unified picture. 

Our results might be of help in different respects. Recently the (T) dualities discussed 
here have played important roles in D-brane and string compactifications and our findings 
concerning the relation between holonomies and constituent positions is an interesting one 
also in this context. For the low charges studied here we saw that the basic constituent is 
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a twisted instanton, with charge Q = ^ in SU{2), which shows up as a BPS monopole at 
finite temperature. We have shown that this fractional instanton maps onto itself under 
the Nahm duality transformation, a property guaranteed by the fact that for SU{2) the 
twist is preserved under the duality transformation and the only freedom is the position of 
these twisted instantons. It would be interesting to see in which sense these constituents 
do continue to play an important role at higher charges. Finally, we hope that the work 
presented here will lead to a more complete analytic understanding. 



Appendix 

In this appendix we will provide the proofs of the basic ingredients of the fiavour construc- 
tion of Nahm's transformation for non-trivial twisted boundary conditions {k, m 7^ 0). In 
particular, we will present in a basis independent way the derivation of the characterisation 
for the Nahm-dual torus and its corresponding dual twist. This complements the results 
of Ref. m. 



Our starting point is the set of Weyl zero- modes \l/^"(a;) satisfying the boundary con- 
ditions of Eq. d^). Consider now the behaviour of these zero-modes under translation of 
the variables z. It is easy to see that given a solution of the Weyl equation one can 

construct a new solution of the Weyl equation as follows: 

;^fe)(a;)^e-2-^(^)vl>,+,(x) , (17) 

for any element y & (as usual the space of linear forms can be identified with R^, and 
y{x) = y ■ x). Applying this to an orthonormal set of solutions of the Weyl equation one 
obtains a new set. This set satisfies different boundary conditions, obtained from Eq. (^ 
by replacing the twist eating matrices T{a), 

r(a) r'(a) = e'"'^(") T{a) . (18) 

We thus see that there is a correspondence between translations in z and different choices 
of twist-eating solutions. 

In general the matrices T'{a) constitute an inequivalent set of solutions to the original 
r(a). However, for special values of y the two sets are unitarily equivalent (related by a 
similarity transformation). For those values of y, denoted by d, there exist U{Nq) matrices 
r(d) satisfying 

f (a) r(a) = e^'^*'^^^) T{a) f (a) . (19) 

This formula shows a nice duality between both sets of U{No) matrices. Before investigating 
for which a these equations have solutions, we want to discuss its consequences. For all 
values d, for which Eq. ( [I9| ) has solutions, we can construct the functions 

e-^-'-(-)xi,il^{x)t'\a) (20) 

which form a new orthonormal set of solutions of the Weyl equation with the same boundary 
conditions Eq. Thus, they can be expressed as a unitary combination of the original 
set of zero-modes. We arrive at the main formula: 

^r+a(^) = e''^^'^") TUa) (x)(fii(z))^" . (21) 
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From this one can deduce that for all a the Nahm transformed gauge field satisfies 

A^{z + a) = [na{z)]A^{z) (22) 

One easily sees that the set of a define a lattice A, which we will call the Nahm-dual lattice. 
It is a sublattice of A, the lattice dual to A. We define the Nahm transformation of the 
original self-dual gauge field with twist as the restriction of A to the Nahm-dual torus 
R^/A. This is a U{NqQ) self-dual gauge field with topological charge Q = NqN/\A/ A\, 
where |A/A| stands for the number of unit cells of A that fit in a unit cell of A, or the 
number of elements in A(mod A). The matrices Cla{z) are the corresponding twist matrices. 

Now let us investigate the solutions to Eq. (|19|) and the characterisation of A. For that 
it is essential to use some properties of the twist eating solutions. One can prove [0] that 
the set of matrices r(a) contains a set of linearly independent matrices, which form a 
basis of the space of all Nq x Nq matrices. Actually, to construct this basis one has to select 
a representative of the quotient lattice A/Aq. Its order is precisely Nq. The sublattice Aq 
can also be defined as the one associated to the matrices r(a) commuting with all other 
r(a') for any a' G A. We call it the center lattice. Using this fact one can, after some effort, 
prove that indeed the matrices r(a) do coincide, up to a phase factor, with the original 
matrices r(a). In more mathematical terms, we can say that there exists a mapping 
from A to A such that 

f(a) = e''^(^)r(0(a)) . (23) 

where (f{a) is an arbitrary phase. From this we can identify A. Let a be an element of Aq, 
then it commutes with all matrices r(a') for any a G A', this includes the matrices r(a). 
Inspecting Eq. ([T9|) we conclude that a (a) G Z for any a G Aq. Therefore A = Aq, the 
Nahm dual lattice coincides with the dual lattice of the center lattice. One can also see 
that the center lattice of A coincides with A. From here it is immediate to see that the 
Nahm-dual of the Nahm-dual gives back the original lattice. 

There is one extra ingredient which has to be worked out, namely what is the Nahm 
dual twist. This follows from Eq. (|23|) and Eq. (H). We have 

f(a) f(a') = e-^-^^^'^'^')) f(a') f(a) (24) 

with 

^(a,a') = -Ar(0(a),0(a')) . (25) 

Eq. ( pTD is only consistent provided the Nahm-dual twist matrices Cla{z) satisfy twisted 
boundary conditions analogous to Eq. (H) with the skew symmetric form A/" replaced by 
Af, the dual twist form. 
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